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Z 2 and Chern topological phases such as newly discovered quantum spin Hall and original quantum 
Hall states hardly both co-exist in a single material due to their contradictory requirement on 
the time-reversal symmetry (TRS). We show that although the TRS is broken in systems with a 
periodically driving ac-field, an effective TRS can still be defined provided the ac-field is linearly 
polarized or certain other conditions are satisfied. The controllable TRS provides us with a route to 
manipulate Z 2 and Chern topological phases in a single material by tuning the polarization of the 
ac-field. To demonstrate the idea, we consider a generic honeycomb lattice model as a benchmark 
system that is relevant to electronic structures of several monolayered materials. Our calculation 
shows that not only the transitions between Z 2 and Chern phases can be induced but also features 
such as the dispersion of the edge states can be controlled. This opens the possibility of manipulating 
various topological phases in a single material and can be a promising approach to engineer some 
new electronic states of matter. 


Introduction. The discovery of topological insulators 
(TIs) in condensed matter systems has not only revealed 
novel physics of the quantum world but also unified 
many physical phenomena, which were thought to be 
irreverent, into the same framework|1. Their peculiar 
edge states make TIs a hot topic for both fundamen¬ 
tal interests and industrial applications. Several materi¬ 
als such as HgTe/CdTe quantum well, Bi^Sbi-^ alloys, 
Bi 2 Se 3 and B^Tes, etc., have been proven to be TIs by 
experiments 2 . Despite these successes, how to design a 
topologically non-trivial material, remains a challenging 
issue. In most cases, the discovery of new TIs still relies 
on serendipity rather than predetermination. 

Instead of searching for materials with intrinsically 
non-trivial topology, there are several recent studies fo¬ 
cusing on manipulation of topological phases using con¬ 
trollable physical processes, e.g. electric fields, strains, 
etc[3H5]. Those studies not only offer new tools to gen¬ 
erate various topological phases but also open new ways 
to making real electronic devices. 

One of the promising methods to engineer a topolog¬ 
ical property of a system is to use periodically driving 
fields 16 fib. The proposal is based on the Floquet the¬ 
ory which states that the Hamiltonian of a system with 
a time-dependent periodic potential can be mapped into 
an effective static Hamiltonian, called the Floquet Hamil¬ 
tonian. If the (quasi)band structure of a Floquet Hamil¬ 
tonian exhibits a topological behavior, we can expect 
there exists a similar feature in the original Hamilto¬ 
nian in a dynamical fashion. An advantage of using this 
method to engineer the band topology is that the ac- 
field provides a set of tunable parameters such that a 
variety of band structures unaccessible in the original 
material can be generated in a dynamical way. Many 
proposals based on the topology of Floquet Hamiltoni¬ 
ans have appeared recently, some of which are: Floquet 
TIs in graphene pT], Floquet TIs in semiconductor quan¬ 


tum wells [12], Floquet Majorana fermions in topological 
superconductors [T3]. merging Floquet Dirac points [14], 
Floquet fractional Chern insulators [15], Floquet Wely 
semimetal [16], etc. A few experiments that support the 
idea of Floquet TIs have also been carried out [17] 118] . 
Those works not only lighten up the road to manipulate 
topological phases but also bring us a vast landscape of 
new physical phenomena that are hardly found in static 
systems. 

While many topological phases have been studied 
within the Floquet framework, the discussion of Z 2 
phases remains scarce because time-reversal symmetry 
(TRS), a necessary condition for the existence of the Z 2 
phase, is always broken due to the time dependence of 
the external perturbation. However, the Floquet Hamil¬ 
tonian is merely an effective mapping of the original 
Hamiltonian, so the loss of TRS in original Hamiltonian 
does not necessarily result in the loss of TRS in Flo¬ 
quet Hamiltonian. Establishing an operator that links 
Floquet states in the Brillouin zone by a similar way as 
conventional TR operator does it for Bloch states, an ef¬ 
fective TRS can be defined!?] fl2] . If so, two seemingly 
contradictory phases, TRS protected Z 2 phase, such as 
recently discovered quantum spin Hall state, and TRS 
broken Chern phase, such as much celebrated original 
quantum Hall state, can both be manipulated in a sin¬ 
gle material by tuning the ac-field, which is the main 
message of the present work. 

Here, we first show how we truncate the Floquet 
Hamiltoian to finite dimension in a realistic calculation. 
Second, we show that the TRS conditions can be eas¬ 
ily satisfied if the field is linearly polarized or certain 
low excitation conditions are reached. Third, we use a 
prototypical 2D material with strong spin-orbit coupling 
as a benchmark in our calculation, in order to demon¬ 
strate the idea of manipulating Z 2 and Chern topological 
phases in the same system. More specifically, we consider 
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a generic half-filled p-orbital honeycomb lattice model to 
illustrate our findings. Our results show the evidence for 
the Z 2 phase with the ability to control the dispersion 
of its edge states by properly tailoring the external ac- 
field. When the polarization is away from the effective 
TRS condition, a rich Chern phase diagram begins to ap¬ 
pear which suggests a Z 2 -Chern phase transition. Thus 
we demonstrate the possibility of manipulating between 
the two topological phases in a single material which can 
serve as a promising tool to engineer some novel elec¬ 
tronic states in condensed matter systems. 

Floquet Theorem. We consider a tight-binding Hamil¬ 
tonian with an external time perodical ac-field H(r) = 

r) + h.c. where r is time, (a, /3) 
are the internal degrees of freedom (e.g. orbitals, 
spins, etc.) of the unit cell positioned at R j and R 
The ac-field is coupled to the problem by introducing 
a minimal coupling t^f 3 (r) —>> ^^ e * A ( r )( r j ~ r f) where 
vj is the position vector of the state | a) in the unit 
cell located at R j and A is the vector potential of 
the field. Since the Hamiltonian has both lattice and 
time translational symmetries, we can perform a dual 
Foruier transform , such as c a j(r) = ^2 n c a j n e~ inujr = 
N~ D / 2 C cm(k)e - 2 ( k '' R j+ nwr ). The Floquet the¬ 

orem proves that the Floquet Hamiltonian Hp{r) = 
H(t) — id r in the Fourier transformed (k, u) space can 
be expressed asp~5l 

Hp{ k,w) = Y, (Ki-nuS nm^ot /?)4n( k )C/3 rri (k) + h.c., 
nmotfi 

hf m ( k) = YKi J m-n(A(r) ■ Ar)]e ik ' R| , (1) 

l 

J g (x(r)) = ^ jF^-^dr, 

where Ar = i*q — rf, k is the wave vector, uj = 2ir/T is 
the frequency of the ac-field and (n, m) are the Floquet 
indexes. 

Because Eq|I] is a key result of the Floquet theorem, 
we assert here that 1 ) Similar to the undriven system, 
the Floquet Hamiltonian Hp forms an eigenvalue prob¬ 
lem Hp(k, cj)|^ 7 n (k)) = e 7 n (k)|^ 7 n (k)) where 7 is the 
band index, n is the Floquet index ranging — oc to +oc 
and e 7n is the so called quasienergy; 2 ) Similar to the 
existence of reciprocal lattice and the periodicity in the 
k-space, the relations e 7n = e 7 o + neo and \u in ) = |^ 7 o) 
are held as a result of the analogous properties of the 
Brillouin zone in the frequency domain. They also show 
the physics of absorbing/emitting n photons, so the Flo¬ 
quet bands are shifted by 3) The solution of the 

original Hamiltonian is obtained by linearly combining 
static Floquet band states |^ 7 (r)) = e~' le ^ r \u 7 (r)) = 
e -ie 1 T ^+^00 e inujT \u in ) where \u 7 (r)) is the Floquet 
state which is periodic both in space and time. Note that 
r no longer appears in Hp and \u in ), so the Floquet the¬ 
orem simplifies the original time-dependent problem by 
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FIG. 1. (color) (a) Floquet bands within first order emis¬ 
sion/absorption photon processes as replicas of the original 
band structure modified by Jo (see left upper inset), (b) for¬ 
mation of the Floquet band structure by merging states into 
single Brillouin zone (left) and accounting for the effect of 
gap opening due to Ji (right), (c) The imaginary part of 
Ji(A • Rj) as a function of polarization (f> (in unit of 7 r) when 
A = [1, 7.467 ]/a (with h = e = 1). Left lower inset gives the 
definition of each R^ of graphene honeycomb lattice, where a 
is the length of R 3 . 


mapping it to a static one. Therefore we can treat Hp 
as the usual lattice Hamiltonian and explore its topology 
using the techniques developed for static systems. If Hp 
has non-trival edge states, we can expect a dynamical 
analogy on |^ 7 (t))[T 2 ]: 4) The form of h^( k) is just the 
usual lattice Fourier transform of the states labeled by 
two indexes (a, n) rather than by one. The extra degree 
of freedom is the penalty of mapping the time-dependent 
Hamiltonian into a static one. The hopping integrals be¬ 
tween the two states c^ n and c^i m are obtained by mod¬ 
ifying the undriven hoppings: t®f —>> t^J 3 J m _ n (A • Ar). 

Because the Floquet index n ranges from — oc to + 00 , 
the Floquet Hamiltonian is not manageable unless we 
make some approximations [ 6 |. Two approximations are 
frequently adopted: (a) weak intensity limit and (b) 
high-frequency limit. 

For the approximation (a), let us consider an ac-field 
sinusoidal in time. In this case, J q ( A • Ar) is essentially 
the g-th Bessel function of the first kind. In the limit of 
the weak intensity, \A\ —>> 0 , its asymptotic behavior is as 
follows: Jo —y 1, Jq^L 0 —» 0. The larger the q the faster 
J q ^0 drops to zero. Hence we can truncate Hp to a finite 
dimension by including just a few lowest order photon 
processes, provided the field intensity is weak enough. 
For example, if we keep q = 0,1, Hp is reduced to the 
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following form 


H f ~ ViHpP- 1 =H) = 

( hf - 1 ul hf 0 \ 

hf hf hf 

\ 0 hf hf + lulj 

where Hj denotes a reduced Floquet Hamiltonian de- 
scribing an emission/absorbtion of a single photon and 
Vi is the operator that projects Hp to Hj . A diagram¬ 
matic explanation of such first order process is shown in 
Fig. 1(a) and (b). In the upper left inset, the undriven 
band structure is modified by the 0 -th order effect Jo. 
Once J\ term comes in, the bands will have three copies 
with energy shifts 0, ±cj. When those bands reach reso¬ 
nant energies, i.e the band crossings, J\ will open gaps 
~ tJ\ making them anti-crossing. This is the main idea 
of the truncation. One has to note that when Hp is trun¬ 
cated, the periodicity in frequency domain is broken, and 
the relations e 7n = e 7 o+ncj and \u in ) = |^ 7 o) are strictly 
speaking no longer valid. 

As for the approximation (b), let us assume the fre¬ 
quency of the external field is so much larger than the 
bandwidth, uj W, that the Floquet bands do not cross 
anymore. In this limit, the gap openings due to J q> o 
become less important, which implies that it is also the 
condition to consider just the lowest order photon pro¬ 
cesses. 

Time-Reversal Symmetry . In an undriven system, the 
TRS is defined by TJ7 (t)T _1 = H(—r) where T is the 
conventional TR operator T = e~ l7T(7y ^ 2 K. Although 
systems with time-dependent ac-fields do not hold this 
property, it is still possible to define an effective TRS for 
the Floquet Hamiltonian[71 fT 2 ). To give specific condi¬ 
tions holding the effective TRS, we conclude with two 
important theorems here (see Supplementary Materials): 

Theorem I: If there exists a parameter To such that 
7~J7 (t)T _1 = H(t + To), one can always define an effec¬ 
tive TR operator Q = c %H fT °T that satisfies the relation 
QHp^Q - 1 =Hp{-k). 

Theorem II: Assuming a system has TRS when it 
is undriven, i.e. |A| = 0 , then A(t) = [A x sin{ujr + 
(j) x ), A y sin((jJT + (j) y ), A z sin(ujT + (j) z )\ with fa — <pj = mir 
(i,j G x,y,z; m G integers ) will automatically make 
H(t) satisfy 7~J7 (t)T _1 = H(r + To). Furthermore if 
the time frame is properly chosen, one can always let all 
fas = nfK (ni G integer) such that To = 0 and Q = XT 

These theorems tell us if the phase differences among 
each field component are multiples of 7 r, the Floquet 
Hamiltonian will have effective TRS [16] and the TR op¬ 
erator can be treated as a conventional one acting in the 
Hilbert space of the basis of the Floquet Hamiltonian 
{|cm(k))}. In the following, we will call the condition 
fa — (f)j = rmr as a linear polarization although the po¬ 
larization is not definable if A is not in 2 D. 



FIG. 2. (color) Honeycomb lattice irradiated by electric ac- 
field A(t). Left inset: band structure without SOC. Right 
inset: band structure with SOC A = 3 1, t being a unit of 
energy. Note that all the bands are doubly degenerate. 


The linear polarization condition is not the only op¬ 
tion to have effective TRS. Since we are handling the 
v-th order reduced Floquet Hamiltonian HJ rather than 
the original Hp in a realistic calculation, it is possible 
that Hj has more time-reversal points than Hp. Recall 
that the hopping integrals in the Floquet Hamiltonian are 
generated by modifying If one 

can properly tailor A such that J q < u are real numbers 
for all Ar in the lattice, obviously TRS will be kept up to 
v -th order HJ. Since Jo is always real, the highest order 
v should be equal or greater than 0. To give an example 
of v > 0, we have plotted in Fig. 1(c) the imaginary part 
of J\ with respect to three non-equivalent position vec¬ 
tors of a honeycomb lattice as a function of polarization 
0 — 4>x — 4>y by setting [A x ,A y ] = [1,7.467]/a. One can 
immediately notice that there are two additional TRS 
points (all lines reach 0) other than 0 = m7r, i.e. 0 = 7t/2 
and <j) = 37t/2. However, those TRS points are just re¬ 
sults of low excitation approximation. One should always 
confirm that the energy splitting A E of Kramer degen¬ 
erate states due to higher order terms is much smaller 
than the characteristic energy e c that we are interested 
in (A E ~ tJ v +1 <C e c ) to explore this feature further. 

Floquet Topological Phases. The best candidates 
to realize the transition between a TRS Floquet 
phase and a non-TRS Floquet Chern phase would 
be 2D materials with spin-orbit coupling, e.g transit- 
metal dichalcogenides| 2 T, graphene with adatoms L IT 
[20] . silicene[ 22 ], germanane[23], Tin films [24], <a-Sn[25]. 
etc. These materials have been proven (or have high 
expectancy) to exhibit monolayer structures with band 
gaps around dozens to hundreds meV. Because of their 
planar geometry, the in-plane ac-field can be easily real¬ 
ized by a laser in experimental setup. 

Here we consider a nearest neighbor tight-binding 
Hamiltonian on a honeycomb lattice with a p-orbital (to¬ 
tal six states) per each site as a generic minimal model 
describing the 2D material at the center of interest. For 
simplicity, we assume the system is half-filled. In order to 
make our model close to actual band structures, hopping 










4 



FIG. 3. (color) Floquet topological phase diagram. Upper 
left: Z 2 phase diagram (green for Z 2 = 0 and blue for Z 2 — 
1) with linearly polarized ac-field [</>x,<l>y\ — [0,0]. Upper 
right: Chern phase diagram (Chern numbers for each color 
are shown in the legend) with eliptically polarized ac-field 
[(j) x , (j)y] = [0.0, 0.5]7r. A x and A y are chosen to be 0.3n/a, n = 
1 ~ 13 (a is lattice constant, h = e = 1). There are four points 
at the phase diagram labeled by 1 to 4 with their edge-state 
band structures shown (all energies in units of t) . Edge states 
are calculated in a zigzag ribbon geometry with 20 sites in the 
transverse direction. 


integrals are generated by a Slater-Koster method 26] 
with V ppa = t, V pp7T = -t and onsite energy E Pxyz = 0. 
SOC is treated as a local potential by evaluating the ma¬ 
trix elements (pi\XL • S| pj) with A = 3 1 for each site. 
In order to calculate the topological invariants, we im¬ 
plement the n-field method introduced Fukui et al. m- 
This method has been proven to provide evaluations of 
both Z 2 and Chern topological invariants in discretized 
Brillouin zones accurately and efficiently. We empha¬ 
size extra time that when computing Z 2 invariants for 
the Floquet Hamiltonian, the TR operator should be re¬ 
placed by the effective TR operator Q = e lHFT °T as 
described in this work. 

In Fig.2, we show a cartoon of the honeycomb lattice 
and the band structures with and without SOC. One can 
find that the Dirac points at high symmetry points T 
and K become gapped when SOC is turned on. This is 
a general feature of a system with the honeycomb lat¬ 


tice. To study the Floquet effect, we consider the re¬ 
duced Floquet Hamiltonian to first order, Hj , and use 
a rather large frequency u = lOt (larger than the band 
width ^ St). The amplitudes of A x and A y are cho¬ 
sen to be 0.3n/a, n = 1 ^ 13 (a is the lattice con¬ 
stant) with linearly polarized field (j) x — = 0- The 

Floquet bands are also assumed to be half-filled as in 
the undriven case. Fig.3 shows the Z 2 phase diagram 
and selected band structures (labeled by 1 to 4) of the 
edge states. Apparently there exists a large (shown in 
blue) area of Z 2 phase in the parameter space. To check 
the corresponding edge states, we have plotted the Flo¬ 
quet band structures of zigzag ribbons under the same 
ac-field. Plots corresponding to the parameters of phase 
points 1 and 3 are both Z 2 = 1 phases, so the Dirac cones 
appear in the gapped region as expected. An interesting 
finding for the phase point 3 is that the Dirac cone ap¬ 
pears in k = 0 rather than at k = n/a as seen for the 
point 1. It is supposed to be the case of the armchair rib¬ 
bon in the Kane-Mele model which now appears in the 
zigzag edge with appropriate parameters [28]. It means 
the ac-field can not only tune the Z 2 topology from triv¬ 
ial to non-trivial, but also control the specific features of 
the edge states. Finally, we have also plotted the band 
structure corresponding to the trivial phase (point 2) as 
a confirmation that the edge Dirac cone indeed does not 
show up. 

Now we discuss how to make the Z 2 phase transiting to 
a Chern phase. Let us consider an elliptically (circularly 
if A x = A y ) polarized ac-field with [(f>x,(/)y\ = [0,7r/2]. 
The phase diagram of the Chern numbers is shown in the 
upper right of Fig.3. Because it is a multiband problem 
(12 bands for the undriven and 36 bands for the reduced 
Floquet Hamiltonian), the Chern number can be much 
larger than ±1[29]. We also show the edge states corre¬ 
sponding to the parameters of the phase point 4. The 
existence of two Dirac cones at the edge agrees with our 
C = — 2 result very well. The phase diagrams of Fig.3 
illustrate how Z 2 and Chern topological phases can be 
manipulated in a single material using properly tailored 
ac-fields. 

Finally, we estimate some physical quantities relevant 
to realization of such exotic electronic phases in real 
systems. Let us take graphene with adatoms as an 
example [T9]. It is predicted to have SOC induced gap 
Eg around 5 ^ 20 meV. To simulate this problem, we 
use tight-binding parameters for the s and p states of 
graphene obtained by fitting to its band structure 30 j 
and tune the SOC to a value that in our model fits 
the gap value of 5 meV. We consider two cases: a mi¬ 
crowave field, uj = 2 GHz <C E g /h, and an infrared field, 
ljJ = ‘ITHz Eg/h. Polarization angles are chosen to 
be 0 and n in order to observe Z 2 and Chern phases re¬ 
spectively. To ensure the weak intensity approximation, 
we limit A x ,A y < 1 (h/A) for both cases so that J 2 ef¬ 
fects are about two orders of magnitude smaller that J\ 
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and can be neglected. The electric field and the corre¬ 
sponding laser intensity are obtained by Eq = Auj/ e and 
Jo — |eoc£o = 1.33 x lO -3 ^ (W/crn 2 ) respectively. For 
the microwave case, we found that the Floquet Z 2 phase 
can be easily observed but the Chern phase cannot. This 
corresponds to the electric fields Eq < 131 Vjcm or the in¬ 
tensities Jo < 23VF/cm 2 , which can be achieved by lasers 
with powers P < 200raLF, easily accessible in experi¬ 
ment. As for the infrared case, we found both Z 2 and 
Chern phases can be generated within that regime. It 
corresponds to the electric fields Eq < 1.31 x 10 5 V j cm 
or the intensities Jo < 2.3 x 10 7 W/cm 2 . This will re¬ 
quire a rather high power about several kW in experi¬ 
ments. This power is still experimentally accessible but 
most materials can burn out under such a strong field. 
Therefore searching for a material that can display both 
phases under lower intensities could be an interesting 
topic for future research. Although realization of Z 2 -to- 
Chern phase transition could be difficult in experiments 
for our discussed system, we have to emphasize that eas¬ 
ily achievable Floquet Z 2 phase still remains a treasure in 
problems of engineering topological electronic structures. 

Conclusion. In summary, we have developed a frame¬ 
work to study TRS in Floquet Hamiltonian and used 
a generic tight-binding model of the honeycomb lattice 
relevant to several recently discovered monolyaered ma¬ 
terials in order to demonstrate the transition between Z 2 
and Chern phases by tuning the polarization of the ac- 
field. Although, our discussion is based on the dynami¬ 
cal analogies, the physics is still very fascinating not only 
due to the emergence of the phase in a formally time- 
reversal breaking potential but also due to the possibility 
of manipulating self-contradictory topological phases in 
a single material. Both phenomena are hard to find in 
static systems but could lead us to a new physics that is 
unreachable in conventional solid-state matter. 
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FLOQUET TIME-REVERSAL SYMMETRY 

Define Floquet operator u(T) and Floquet Hamilto¬ 
nian Hf 

u(T ) = T[e _i J° Tff ( r)<iT ] = e~ iHpT 

where T is the time periodicity of the Floquet system and 
T is the time-order product. We hope to find an effective 
time-reversal (TR) operator Q for the Flouqet operator 
and Floquet Hamiltonian such that 

QHFQ- 1 = u(-T) 

and 


origin from ro to 0, then an effective TR operator can be 
defined as Q = 1ZJC 

1Zu* (T, 0)7^ = u( 0, r 0 )Su* (r, 0 )<S^t (0, r 0 ) 

= w(0,r 0 )w t (r 0 + T,To)w t (0,r 0 ) 

= u\T,0) 

It means 

Qu(T,0)Q~ 1 = u(—T,0) 

and 

QH F (\a)Q~ l = H f (- k) 


QHp^Q - 1 = H F (-k) 

where Q is an antilinear oparator with Q 2 = —1. We 
claim that if there exists an parameter To that satisfies 
the relation 

TH{t)1 1 = H(—t + T 0 ) 

where T is the conventional TR operator. Then an ef¬ 
fective Q can always be defined as 

Q = u(0,T 0 )T=e iHFTo T 

In the following, we provide a proof for this theorem. 
(Note: Our proof is equivalent to the one shown in Ref.7 
of the main article. Because we have chosen a slightly 
different statement, we prove it again here.) 

Let us represent the conventional TR operator as the 
product of an unitary operator S (usually e ~ l7T(7y / 2 ) and 
the complex conjugate operator JC: 

T = SJC 

Assume there exists a parameter To such that 
Ti?(r)T _1 = H(-t + t 0 ) «=>■ SH*{t )<S f = H\-t + t 0 ) 
Since 

u(T, 0) = lim e -' lArH ( T ~ Ar ) x e -^ rH ( T - 2Ar ) x ... 

N^oo 

• • • x e“* ATff (0) ; At = T/N 

We have 


RELATION TO POLARIZATION 

Consider a vector potential A(t) = [A x sin(ujr + 
(j) x ), A y sin (cor A <fi y ) , A z sin (ujt + <fi z ) ]. We claim two con¬ 
sequences: 

• If <fii — (fij = mir (■ i,j G x,y,z, m G integer ), 
the Floquet time-reversal criterion: T77 (t)T _1 = 
H(—t + To) will always be satisfied. 

• If <fii — (fij = mir, one can always let <fii = n^ (n^ G 
integer) such that To = 0 and the effective TR 
operator can be simply expressed as Q— TT ’. 

The first theorem tells us the relation between the Flo¬ 
quet TR symmetry and the polarization of the ac-field. 
The second theorem helps us to deal with the effective 
TR operator in a much simpler way. In the following, we 
provide a proof for these two statements. 

Consider the basis set of Hilbert space {|cqcr)} where 
a is the label of space-related degree of freedom, e.g, 
sublattice, orbital, etc., and a = +/— is the spin index. 
Define time-reversal operator T = uK where u = — ia y . 
Then the matrix element of a time-reversal transforma¬ 
tion applied to the Hamiltonian is given by 

(acr\THl 1 |<a a ) = (a<j\uH*u^\a a ) 

= (-1 \l 5 <r-+ 5 a '\* 

' ' ' a —a ' 

If TR symmetry exists, THT~ l = H , and we obtain a 
restriction on the matrix elements: 


Su'fr.OW* = lim e i^rSH-(T-Ar)S^ x ... x ArSH* (0)5+ 

N^oo 

— \[ m (j jAr (to —T +At) x . . . x giAr iT (r 0 ) 

AT—>• oo 

= lim e ^r^ t (r 0 +Ar) x . . . x g iA tH\t 0 +T) 
iV—>- oo 

«'(-(, i-Y'.r (1 ) 

where the third equal sign uses the relation H(r + T) = 
H(r). Therefore if we define 7 Z = u(0,tq)S to shift the 


' ' ' rv' — rr' ' CM. (7 


For a system with an ac-field, the hopping integral is 
modified by t a 7 ,(r) t a 7 ,e zA( ' r ^ rQi_r « / \ The Floquet 

J a a v J a a 1 

TR criterion TH{r)T~ l = H(—r + tq) requires that 


(_l )[ <5 ^- +(5 a'-\(^t a ~ a / )*e _iA( ^ r ^ ra_r a 


= t a 7 / e iA ( _r+r °^ ra_ V) 

a a 
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Assuming the system has TR symmetry when un¬ 
driven, the hopping integrals will be canceled out and 
we have 


- A(t) = A(— T + To) (2) 


Because 


A(r) = [AxSin^ojr-^-fpx), A y sin{ijjT-\-(l)y), A z sin(ujT-{-(/) z )\ 


then Eqj2] means 

— (UJT + (j)i) + 2n^7T = —ujt + UJT o + <fi x (3) 

=>* 4>i = riiTT — cjto/ 2, i G x,y,z ; rii G integer 

Since To can be arbitrary real numbers, it is convenient 
to express the effective TR condition as 

ct>i - (j)j = m7r, rri E integer (4) 

For an in-plane ac-field, it is to say that a linearly polar¬ 
ized ac-field will have the Floquet TR symmetry. Fur¬ 
thermore, if Eqj4] is held, we can always choose To = 0 
and Q = XT . To show this, let us shift the time frame 
t = t' Tto/ 2 and plug in the Eqg If so, we can get a 
new equation where <pi = ni 7r and t 0 = 0. 


